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Abstract 

e-j \ We have analyzed 200 Rossi X-ray Timing Explorer observations of the black hole candidate GX 339-4, 

O |. all from the bright hard state periods between 1996 and 2005. Purpose of our study is to investigate the 

radiation mechanisms in the hard state of GX 339-4. The broadband 3-200 keV spectra were successfully 
modeled by a simple analytic model, power-law with an exponential cut-off modified with a smeared edge. 

■ The obtained energy cut-off (E cut ) was distributed over 50-200 keV, and the photon index over 1.4-1.7. 
We found a clear anti-correlation (i^cut L _0 - 70±0 06 ) between the X-ray luminosity (L) in 2-200 keV and 
E cut , when L is larger than 7 x 10 37 erg s _1 (assuming a distance of 8 kpc), while £ cut is roughly constant 
at around 200 keV when L is smaller than 7 x 10 37 erg s _1 . This anti-correlation remained unchanged by 

. adopting a more physical thermal Comptonization model, which resulted in the anti-correlation that can 

' be expressed as kT e oc £-°- 24 ± 06 . These anti-correlations can be quantitatively explained by a picture in 

0^ ■ which the energy-flow rate from protons to electrons balances with the inverse Compton cooling. 

r*** ' Key words: accretion disks, black hole physics — stars: individual (GX 339-4) — X-rays: stars 

■ 1. Introduction fully modeled by thermal Comptonization model (Dove et 
' al. 1997, Pountanen & Svensson 1996), although location 

Black hole X-ray binaries show large varieties in their and geometry of the Comptonizing medium is still de- 
X-ray properties. Several spectral and variability states bated. Other models include advcction dominated accre- 
have been recognized in the past, with varying names and tion flows (ADAF) (Esin et al. 1997), synchrotron and/or 

■ varying definitions (see e.g. McClintock & Remillard 2006, synchrotron self-Compton radiation from the base of a jet 
?H \ Homan & Bclloni 2005) and it is believed that they cor- (Markoff, Nowak, & Wilms 2005), and external Compton 

respond to different accretion geometries. In this paper scattering in the jets (Georganopoulos, Aharonian, and 

we focus on the spectral properties of the so-called hard Kirk 2002). 

state, which is the state observed at low X-ray luminosi- GX 339-4 was discovered with the X-ray satellite OSO- 

ties (i.e. below a few percent of the Eddington luminos- 7 in 1971 (Markert et al. 1973). Since it was similar to 

ity, £Edd), although it can also be seen during the ris- Cygnus X-l in terms of X-ray spectral and variability 

ing phase of transient outbursts below ~0.2iEdd- The properties, it was considered a black hole candidate (BHC, 

accretion geometry in some of the X-ray states is rel- Samimi et al. 1979). A recent measurement of the mass 

atively well understood - e.g. the soft state is thought function (5.8±0.5 M Q , Hynes et al. 2004) strengthens 

to be governed by an optically thick and geometrically this. We thus assume a mass of 5.8 M Q for this black 

thin accretion disk - this is not the case for the hard hole candidate as a secure lower limit. GX 339-4 is one 

state. X-ray spectra in the hard state are well repre- of the best-studied BHCs in X-rays and gamma-rays; it 

sented by a power-law with a photon index of 1.4~1.7 was observed with Ginga/LAC (Ueda, Ebisawa, and Done 

(Tanaka & Shibazaki 1996) and, in some sources, a high 1994), CGRO/OSSE (Grabelsky et al. 1995, Smith et al. 

energy cut-off at -100 keV (Grove ct al. 1998). Several 1999), ASCA (Wilms et al. 1998), RXTE (Smith et al. 

models have been proposed to explain the hard state spec- 1999, Belloni et al. 2005), Bcppo-SAX (Corongiu et al. 

tra. Broadband X-ray energy spectra of a large number 2003) and INTEGRAL (Bclloni et al. 2006, Joinct et al. 

of black hole binaries in the hard state have been success- 2006). GX 339-4 has shown various spectral states during 
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Fig. 1. RXTE/ASM light curves (upper panels) and hardness ratios HR2 (= 5—12 keV/3-5 keV; lower panels) of the black hole 
candidate GX 339-4 in the active periods from 1996 to 2005. The arrows labelled a to h indicate the hard state intervals of which 
the RXTE/PCA observations were analyzed in this paper. 



its prc-RXTE outbursts (Tanaka & Shibazaki 1996), but 
often stayed in the low- luminosity hard state (low/hard 
state). All the spectra in the low/hard state can be 
roughly described by a Comptonization model, requir- 
ing the addition of a reflection component, soft excess 
and iron-K lines. Detailed broadband analysis in the 
2-1000 keV range was made using Ginga, RXTE and 
OSSE data from 1991 and 1996 observations by Zdziarski 
ct al. (1998) and Wardzinski ct al. (2002). Wardzinski 
et al (2002) concluded that the four hard state spec- 
tra have very similar intrinsic photon indices of = 1.75. 
Furthermore, they found that the energy cut-off possi- 
bly decreases with increasing the luminosity. Zdziarski 
et al. (2004) studied the long-term behavior by compil- 
ing GINGA/ ASM, CGRO/BATSE and RXTE/ASM data 
(spanning 16 years) and also found that the electron tem- 
perature depends on the luminosity by showing a positive 
correlation between the BATSE flux and the photon in- 
dex in the 70-160 keV range which reflects the breaking 
energy. 

GX 339-4 has also been extensively studied at other 
wavelengths. The optical counter part was first identified 
by Doxsey et al. (1979) as a V~18 blue star. However, 
even during the X-ray off state, optical spectroscopic ob- 
servations with the VLT revealed no spectral features from 
the companion star (Shahbaz et al. 2001). Hence, the 
type and distance of the secondary star are still unknown. 
There are two papers which give a lower limit on the dis- 
tance: Maccarone (2003), 7.6 kpc, and Hynes et al. (2004), 
5.6 kpc. A careful study of the distance is presented by 
Zdziarski et al. (2004) who argue for a most likely distance 
of 8 kpc. We assume a distance of 8 kpc in this paper. 
Multi-wavelength observations have shown that the ra- 
dio emission of GX 339-4 correlates very tightly with the 
X-ray fluxes over more than two orders of magnitude, sug- 
gesting that the jet plasma may also play a role in high en- 
ergy band (Corbel et al . 2000, Corbel et al. 2003, Markoff 
et al. 2003). Corbel et al. (2003) showed that a significant 
fraction of the X-ray flux observed in the low/hard state 
of black hole candidates may be due to optically thin syn- 
chrotron emission from the compact jet. Also jet emission 
have been found in the near-infrared (Corbel & Fender 



2002) and Homan et al. (2005) discovered a tight relation 
between near-infrared and X-ray fluxes in the hard state of 
GX 339-4, similar to that of the X-ray/radio correlation. 

Many authors have only studied a small number of ob- 
servations during individual outbursts of GX 339-4. In 
order to improve our understanding of the radiation mech- 
anisms in the hard state of GX 339-4, we have performed 
a systematic study of detailed correlations among spec- 
tral parameters, using the large archive of public RXTE 
data for GX 339-4. In Section 2 we describe our data 
analysis, and in Section 3 we present our results, focusing 
on the correlations between spectral parameters such as 
luminosity and high energy cut-off or electron tempera- 
ture. Finally in Section 4, we will discuss the origin of the 
observed correlations. 

2. Observation and Data Reduction 

GX 339-4 has been observed for over 10 years with 
the Rossi X-ray Timing Explorer (RXTE, Bradt et al. 
1993). RXTE carries three scientific instruments: the 
Proportional Counter Array (PCA: Jahoda et al. 1996), 
the High Energy X-ray Timing Experiment (HEXTE; 
Rothschild et al. 1998), and the All-Sky Monitor (ASM: 
Levine et al. 1996). 

Figure 1 shows the 1.5-12 keV RXTE/ASM light curves 
of GX 339-4 from three active periods in the last 10 years: 
in 1996-1999, 2002-2003, and 2004-2005 (data from the 
2007 outburst were not analyzed). Following McClintock 
& Rcmillard (2006), we also plot the ASM hardness of 
HR2 (defined as the ratio of ASM counts in 5-12 keV to 
counts in 3-5 keV), where values of ~1.5 are typical in 
the hard state. 

Selection of the RXTE/PCA observations for our anal- 
ysis was based on the location of these observations in a 
hardness-intensity diagram, as shown in Figure 2. This 
diagram was created from the background subtracted 
Standard 2 mode PCA data (see below), where count rates 
are taken from the 3.5-20 keV band and the hardness de- 
fined as the ratio of count rates in the 10-20 keV and 
3.5-5.5 keV bands. All count rates were normalized to 
those of the Crab, to account for possible changes in the 
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detector response. Hard states correspond to the slightly 
curved vertical branches in the right part of the diagram. 
For selecting hard state observations for our analysis we 
used two criteria: hardness >1.4 (to select hard states) 
and count rate > 15 mCrab (to get enough counts per 
spectrum). Thus, we obtained 200 RXTE observations. 
The periods during which these observations were made 
are indicated by the horizontal arrows in Figure 1. A list 
of all the selected observation IDs is given in Tables 2-5. 




0.5 1 1.5 2 

Hardness 

Fig. 2. Hardness-Intensity diagram of the 1996/1997, 
2002/2003 and 2004/2005 outbursts, created from 
RXTE/PCA data. The area in the upper right (outlined 
with the dashed lines) indicates the selected observations. 

For the spectral data reduction we used the HEADAS 
6.0.2 software package provided by NASA/GSFC. 
Additional selection criteria for both PCA and HEXTE 
are as follows: 1) elevation angle from the Earth was 
larger than 10 degree, and 2) the offset angle was smaller 
than 0.02 degree. PCA spectra were only extracted from 
Proportional Counter Unit 2 (PCU2), since it was the only 
PCU that was active in all the selected observations. We 
used Standard 2 mode data, which has a time resolution 
of 16 sec and 129 energy channels covering the full PCA 
(~2-60 keV) energy range. The spectra were deadtime 
corrected (a few %), and the background, estimated by 
the background model for bright sources, was subtracted. 
As for HEXTE, we used the archive mode data with a 
time resolution of 16 sec taken from both cluster A and 
B, and subtracted the background taken from the rocking 
motion. Deadtime correction for extracted spectra was 
made using hxtdead. 

We added 2% systematic errors to each PCA spec- 
tral bin so as to obtain x„ 2 ~ 1 for the spectrum of the 
Crab. The PCA energy range was limited to 3-20 keV and 
the HEXTE energy range to 18-200 keV. Our PCA and 
HEXTE response files gave consistent values of the pho- 



ton index of the Crab spectrum for the PCA (2.08±0.02) 
and HEXTE (2.10±0.02), indicating that the PCA and 
HEXTE could be fitted simultaneously. We did so by al- 
lowing the normalization of the HEXTE spectrum relative 
to that of the PCA spectrum to vary. We found that the 
normalization factor varied between 0.90-0.95. 

Two simple spectral models were applied to all the 
data in carrying out the systematic spectral analysis 
of the continuum. One is a simple analytic model, a 
power-law with an exponential cut-off (cut-off power-law 
in XSPEC, hereafter CPL), and the other is a thermal 
Comptonization model with a spherical geometry, intro- 
duced by Sunyaev & Titarchuk (1980) (compst in XSPEC, 
hereafter COMPST). In addition to their normalization 
parameters, both models are described by two other pa- 
rameters: photon index (a) and energy of spectrum cut- 
off (i?cut) hi the former model, and electron temperature 
(kT c ) and Thomson optical depth (t) of a high tempera- 
ture plasma in the latter. These continuum models were 
further modified by galactic absorption (wabs in XSPEC) 
and a smeared edge (Ebisawa et al. 1994). The hydrogen 
column density in the wabs model was fixed to 5 x 10 21 
enr 2 (Ilovaisky et al. 1986). Edge energy and width in the 
smeared edge model were fixed at 7.11 keV (corresponding 
to the neutral iron-K edge) and 10 keV, respectively. 

Absorption-corrected X-ray luminosities (L) were cal- 
culated in the 2-200 keV range and were based on the 
PCA fit parameters. We assumed the distance of GX 
339-4 to be 8 kpc (Zdziarski et al. 2004). The obtained 
luminosities ranged from l.OxlO 37 erg s^ 1 to 2.1 xlO 38 
erg s _1 , corresponding to 1.3-29% of the Eddignton limit 
(£Edd=7.3xl0 38 erg s -1 ), for a black hole mass of 5.8M Q . 

3. Analysis and Results 

To show the presence of the high energy cut-off in GX 
339-4 in a model independent manner, we calculated the 
PHA ratio of GX 339-4 to the Crab Nebula for several 
observations. A Crab spectrum from 2002 April 28 was 
extracted in the same way as described in the previous 
Section. This spectrum could be fit in the 3-200 keV 
range by a simple featureless power-law model with a pho- 
ton index of 2.10±0.01. GX 339-4 spectra were taken from 
observations on April 3rd 2002 (MJD 52367), April 10th 
2002 (MJD 52374), and April 26th 2002 (MJD 52390). 
Figure 3 shows the ratio of the three GX 339-4 spectra to 
the Crab spectrum. As can be seen from this figure, there 
is a clear cut-off at high energies that changed from ~30 
keV to ~70 keV with decreasing flux, while the spectral 
slope, i.e. photon index, does not seem to vary signifi- 
cantly. 

Next, we performed a systematic study of the spectra 
with the two models described in Section 2. In the CPL 
model fitting, we set the maximum value at 500 keV for 
the energy cut-off -B^ut- The average and standard devi- 
ations of the reduced chi-squarcs (x„ 2 ± a x 2, for 74 de- 
grees of freedom) were 1.0±0.2 and 1.0±0.3 for the CPL 
and COMPST models, respectively. Figure 4 shows exam- 
ple fits to broadband spectra of selected five observations. 
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Fig. 4. Spectral fit results of the combined PCA and HEXTE spectra of GX 339—4 for 5 representative observations. The top panels 
show an unfolded vF v spectrum fitted with a cut-off power-law model, while the three bottom panels show residuals with a simple 
power law, cutoff power-law, and thermal Comptonization (COMPST) model, respectively. 
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to 1.6), while other epochs show no significant variations 
around 1.6-1.7). 
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Fig. 3. PHA ratio of GX 339-4 to Crab taken during differ- 
ent epochs. The high energy cut-off depends on the X-ray lu- 
minosity, while the slope of this spectrum remains unchanged. 

Both models could fit the 3-200 keV spectra for each ob- 
servation reasonably well. The effects of a possible soft 
excess and the iron-line features at 6.4 keV were negligi- 
ble and Xf 2 did not improve significantly when we added a 
disk blackbody or a Gaussian component to our fit model. 
However, in cases when kT e was smaller than 20 keV, we 
found that it was difficult to fit the spectrum with the 
COMPST model, resulting in xJ 2 values that exceeded 2 
(as we show later, reflection models were required in those 
cases). Fit parameters in 5 representative observations are 
given in Table 1, and all the fit parameters are given in 
Tables 2-5. 

The top panels of Figures 5 and 6 show the relation 
between some of the spectral parameters obtained with 
the CPL model. The different plot symbols correspond 
to intervals a—h as defined in Section 2 and as shown in 
Figure 1. The high energy cut-off (Scut) ranged from 40 
keV to 200 keV or more. Although we could not constrain 
the cut-off well at the high energy end, the highest val- 
ues we measured are still well below the upper bound of 
500 keV. As one can see from Figure 6, we found a clear 
anti-correlation between L and S cut when L > 7 x 10 37 
erg s _1 . This relation follows the approximate relation 
Scut oc Z,-0-70±o.(B On the other hand, S cut seemed to 
be roughly constant at 200 keV when L < 7 x 10 37 erg 
s _1 . To confirm whether this energy cut-off was really 
required in the HEXTE data of E~200 keV, we also com- 
pared with the simple power-law fit results. The F values 
(=(Ax 2 /Az/)/(x 2 M)) arc 23.8, 67.4, and 183 for {Au,u 2 ) 
of the upper three observations shown in Table 1, which in- 
dicates that high energy cutoff is required at a confidence 
level of > 99.9 %. The periods when the X-ray luminosity 
was larger than 7xl0 37 erg s" 1 correspond to the initial 
rising phase of the outbursts in 2002 (epoch d) . Unlike the 
high energy cut-off, the photon index was not dependent 
on the X-ray luminosity. It was distributed over 1.4-1.7, 
which is a typical value for black hole candidates in the 
hard state. In the rising phase of the outbursts (epoch d 
and f-g), the photon index becomes steeper with time (1.4 
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Fig. 5. Top: Relation between the luminosity in 2-200 
keV(L) and the photon index(o) obtained by the cut-off pow- 
er-law model. Bottom: relation between the luminosity(L) 
and Thomson optical dcpth(r) obtained by COMPST model. 
The different symbols correspond to the epoch defined in 
Figure 1. 

In the bottom panels of Figures 5 and 6, we further show 
the correlations between some of the parameters obtained 
from fits with the COMPST model. The electron tem- 
perature kT a and Thomson optical depth r were between 
~20-30 keV and ~4-5, respectively. We found a clear 
anti-correlation between kT e and L, kT c oc £~°- 24±0 - 06 , 
when L > 7 x 10 37 erg s _1 , the same luminosity as was 
found with the CPL model. The electron temperature was 
more or less constant around 26 keV when L < 7 x 10 37 
erg s . No clear luminosity dependence was found for r. 

So far we have fitted the spectra with relatively simple 
models. In other works Compton reflection components 
have often been taken into account as well. Therefore, we 
also checked for a number of our observations how the fit 
parameters were affected by using a model for reflection 
from a cold disk (reflect in XSPEC, Magdziarz & Zdziarski 
1995) plus iron-K emission line with a narrow line width of 
50 eV, instead of smeared edge models. The high energy 
cut-off (S C ut) and photon index (a) became systematically 
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Table 1. Best-fit parameters for fits with various model. 



Observation ID 


MUUiiL 


L{w erg s ) 


i/cut / ki c (keV) 


a / r 


Si/Z7T 


X /a.o.i 


i a i a a a i ao a a 

40108-01-03-00 


power-law 


2.2+0.1 




1 TA 1 A AA 

1.70+0.02 




82/80 




CPL 


2.0+0.1 


a r a + 144 

250I 69 


-1 AA 1 A A i 

1.62+0.04 




60/79 




CPL+refl 


2.0+0.1 


>500 


1 C> O I A AO 

1.83+0.03 


n qq+0.12 

U-^-O.lO 


76/78 




COMPbT 


1.9+0.1 


27+3 


j -, o+0. 26 
4 - 18 -0.28 




52/79 




COMPST+rcfl 


1.9+0.1 


.» +29 

46_ii 


__+() 5(, 

z - '°-1.07 


0.26+0.13 


72/78 




COMPPS 


1.8+0.3 




r.oo_ 35 


19+ 011 


87/78 


80132-01-10-00 


power-law 


3.1+0.1 




1.54+0.01 




143/75 




CPL 


3.1+0.2 


A . ) J —1— 1 

2341^ 


1 1 I A AO 

1.41+0.03 




75/74 




CPL+refl 


3.1+0.1 


>500 


1^1 +0 02 

1 -61 -0.04 


n oo+O 03 

0.33ZoIio 


79/73 




COMPST 


2.8+0.1 


30+2 


■ „_l_n 1 7 

4.46l8;12 




98/74 




COMPST+refl 


3.0+0.1 




3.601°;^ 


0.28+0.10 


59/73 




COMPPS 


3.0+0.2 


66+2 


<2.42 


A 1 4+u.uo 


69/73 


80132-01-17-01 


power-law 


6.6+0.1 




1.62+0.01 


249/75 




CPL 




I68I20 


1.44+0.02 




72/74 




CPL+refl 


6.7+0.2 


>500 


w^4-0 02 
' u -0.03 


u -* z -0.04 


104/73 




COMPST 


6.0+0.2 


27+2 


4.57+0.14 




90/74 




COMPST+refl 


6.2+1.1 


36t^ 


+U.28 

°' ul -0.30 


O7+O. 04 

u * z ' -0.08 


60/73 




COMPPS 


6.2+0.3 


58Ii 


<2.49 


„ . ,--,-(-0 0^ 


86/73 


70109-01-02-00 


power-law 


8.1+0.1 




1.64+0.01 




862/74 




CPL 


o i+0 3 


1 1 r+8 
H5_7 


1 O/"* 1 A AA 

1.36±0.02 




80/73 




CPL+refl 


8.4+0.3 


>470 


1 n a 4-0 04 


46+ 09 

u -^°-0.08 


201/72 




COMPST 


7.4+0.2 


25+1 


4.76+0.09 




141/73 




COMPST+refl 


7.4+0.3 


25+1 


. ryt+O 10 

4- 74 l .09 


A OO 1 A A >~7 

0.33+0.07 


1 Or7 / WA 

127/72 




(jUMrrs 


q 1 _i_n 1 
o.l±U.l 


t;a+ 2 
00 _ 1 


<z.yz 


n qq+0.05 
u -°°-0.04 


1 nn /70 
1UU/ ( z 


70110-01-02-00 


power-law 


17.8+0.3 




1 76+ 01 
i - '"-0.02 




534/74 




CPL 


16.8+0.8 


80l^ 


1.44+0.03 




82/73 




CPL+refl 


16.8+0.7 


ziz -40 


1.72+0.05 


41+° 13 


115/72 




COMPST 


16.1+0.8 


20+1 


5.02+0.16 


126/73 




COMPST+refl 


16.3+0.4 


25±| 


4 19+0.30 


31+ - 13 

u - ol -0.12 


81/72 




COMPPS 


16.7+0.5 


47+J 


<2.81 


0.44+0.07 


86/72 



Notes. Errors are quoted at statistical 90% level. 



higher and steeper, respectively, compared to their origi- 
nal fit values (sec Tabic 1 and top panel of Figure 7). For 
example, we derived 1.72+0.05 (original value: 1.44+0.03) 
and 212I4J (79.9lg y) as the new values of photon index 
and E cu t, respectively for observation ID 70110-01-02-00. 
We could not constrain E cut in the low luminosity region 
(L < 7 x 10 37 ) erg s _1 because the reflection component 
strongly couples with a high energy cutoff. The same pro- 
cedure was also applied to the COMPST model fitting. 
The reflection fits gave better results than COMPST with 
a smeared edge model in some observations, when kT e was 
lower than 20 keV, although the electron temperature be- 
came slightly higher. Figure 7 shows the luminosity and 
cutoff/electron temperature relation in the epoch d and g, 
which showed a clear anti-correlation in the non-reflection 

— 1 98 +0 ™ 

fits. The anti-correlation was still found E cut oc L~ l vo -° ™ 
or kT e oc £-°- 18 ± - 05 . These correlations have different 
slopes from other two models, but change still seems to 
occur around the same luminosity, even when we applied 
the reflection component, implying that the change in be- 
havior above L = 7 x 10 37 erg s _1 is model independent. 
The corresponding electron temperature distribution 



in the non-reflection fits is in a narrower range, 16-40 
keV. These values should be taken with some care be- 
cause COMPST model gives a smaller electron temper- 
ature than other Comptonization models, such as the 
COMPPS model (Poutanen & Svensson 1996). We have 
used the COMPST model, the simplest and computation- 
ally fastest of the various Comptonization models, but we 
also tried to apply the COMPPS model with a spheri- 
cal geometry (more complex and accurate) to part of our 
data sets for the five observations covering almost all of 
the luminosity range from l.OxlO 37 to 2.1xl0 37 erg s _1 . 
Table 1 shows the best-fit parameters for them. The range 
in electron temperature that we found was ^40-120 keV, 
which is consistent with the 50-100 keV range found by 
Zdziarski ct al. (2004). In addition, the electron temper- 
ature in the bright hard state [L >7x 10 37 erg s _1 ) is 
significantly lower than that in the low hard state. 
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Fig. 6. Relation between luminosity (L) in 2-200 keV and 
energy cut-off (i? C ut) obtained by CPL model or electron 
temperature(fcT e ) by COMPST model. The dotted lines show 
the best-fit relations for L — E cu t and L — kT e using data in 
the range of L >7xl0 37 erg s" 1 . 



Fig. 7. Relation between luminosity in 2-200 keV(L) and 
high energy cutoff/electron temperature after adding a reflec- 
tion component to out cut-off power-law model. The dotted 
lines show the best fit relations for L — E cu t and L — kT e using 
the data in the range of L >7xl0 37 erg s . 



4. Discussion 

4-.1. Summary of Results and Comparison with Previous 
Results 

We have analyzed a large number of RXTE spectra of 
the black hole candidate GX 339-4 when it was in a bright 
hard state, and have found a clear luminosity dependence 
of the high energy cut-off in the X-ray spectrum. The 
2-200 keV X-ray luminosities rang ed from 1.0 x 10 37 to 
2.1 x 10 38 erg s _1 , covering a factor of 20. This cor- 
responds to 1.3-29% of Eddington luminosity for a 5.8 
solar mass object (assuming a distance of 8 kpc). This is 
the first systematic and quantitive spectral study over a 
wide-luminosity range in the hard state of GX 339-4. We 
further notice that the anti-correlation between X-ray lu- 
minosity and energy cut-off was only found above a lumi- 
nosity of 7xl0 37 erg s _1 . In GX 339-4, such luminosities 
were seen mainly during the rising phase of the 2002/2003, 
when the source was in a bright hard state. This state was 
difficult to observe prior to RXTE, because typical black 
hole transients pass through this state within a few weeks 



of the start of an outburst. RXTE's fast response and 
flexible scheduling were critical in observing this state in 
an increasing number of transient BHCs. 

A high energy cut-off around 100 keV has been observed 
in the energy spectra of many BHCs in their hard state 
(Grove et al. 1998), although it is not observed in the hard 
state of every source as recent INTEGRAL observations of 
GRO J1655-40 show (Caballero-Garcia et al. 2007). For 
GX 339-4 we reconfirmed that it had a high energy cut-off 
in almost all of its hard state observations. It is difficult to 
make comparisons with previous results (Zdziarski et al. 
1998), because the authors used different models, but we 
clearly found a varying energy cut-off in GX 339-4, with 
values between 40 keV and 200 keV or more. High energy 
cut-off energies above 200 keV are beyond the HEXTE 
bandpass, so the highest cut-off energies that we found 
may have some intrinsic uncertainties since a part of them 
were not constrained. Other high sensitive observations 
with INTEGRAL or Suzaku are needed to confirm cut-off 
energies above 200 keV in GX 339-4. 
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4-2. Anti- correlation between Luminosity and 
Temperature 

We found that the high energy cut-off or electron tem- 
perature strongly depended on the X-ray luminosity when 
the luminosity was higher than 7xl0 37 erg s _1 , i.e. O.f 
^Edd- A possible anti-correlation in GX 339-4 was al- 
ready suggested by Wardzmski et al. (2002) and Zdziarski 
et al. (2004). A simple interpretation of the observed 
anti-correlation, i.e. that the hot electrons are efficiently 
cooled by soft photons via inverse-Compton scattering, is 
discussed in more detail in the Section 4.4. Similar anti- 
correlations have already been suggested in GS 2023+238 
(Inoue 1994), GRO J0422+32 (Kurfess 1996) and 10 other 
BHCs including XTE J1550-564, GS 1354-645, and GX 
339-4 (Yamaoka et al. 2005). Hence, this anti-correlation 
might be not intrinsic to GX 339-4 but universal among 
BHCs. When the GX 339-4 luminosity was lower than 
7xl0 37 erg s" 1 , the source seemed to show different be- 
havior; the high energy cut-off and electron temperature 
might have a very weak dependence on the luminosity, 
remaining a nearly constant at 200 keV and 26 keV, re- 
spectively. The reason why the high energy cutoff reaches 
constant value is beyond the scope of this paper, but some 
other cooling mechanisms such as bremsstrahlungs and 
synchrotron emissions may have to be taken into account. 
The hard state of Cygnus X-l, which has a luminosity 
of 0.01-0.03 LEdd, belongs to the latter constant-cutoff 
regime. This may be the reason why this source has 
almost the same cut-off in the hard state at any time 
(Gierlinski et al. 1997). 

4-3. Cooling Time Scale due to Inverse Compton 
Scattering 

Based on the spectral fit results, the thermal 
Comptonization model gives a reasonably acceptable fit. 
In the following, we show whether the Inverse Compton 
scattering can be really effective emission process using 
derived physical parameter. We consider a spherical hot 
plasma (like an ADAF for example) with a radius of 
R (~ 10i? s = 171(g-|^) km; R s : Schwarzchild Radius) 
around a black hole. In this situation, protons will be 
heated up by the gravitational energy release of the ac- 
cretion flow through a mechanism such as viscous heat- 
ing. The protons in the plasma lose their energy due to 
Coulomb collisions with electrons, while the heated elec- 
trons cool due to the Inverse Compton scattering. 

We define two time scales: t pc (oc " lp ra " lc — ; n is 
number density of the plasma, kT p is proton temperature 
and TO p is the proton mass.) as the time to establish a 
Maxwcllian velocity distributions for protons interacting 
with electron field particles (for see details Spitzer 1962), 



given by 



and t, 



comp (- 



-) as a cooling time scale of the elec- 



tron due to inverse Compton scattering. The cooling rate 
per unit volume (^) comp is approximately given as 
E^-ad?"! o"t where c is the light speed, C/ ra d is the photon flux 
density, which is given by U ta ^ — > and <7t is the cross 
section for Thomson scattering. These two time scales are 



"pe 



2.5 x io-4 l3Q fce y MW [sec] 



and 



5.7 x 10" 



(i) 



r_R_y 



( 



10 37 



=*)(§)' 



sec 



(1) 



(2) 



respectively. On the other hand, i v i s (~ ~~§sn )» ^ nc v i scous 
time scale in the ADAF, that is, the time scale in which 
electrons fall to the central black hole, is 



'2.6x 10 



(R\ 2 ( R \ 3 / 2 



(3) 



where 2H is the thickness of disk, fl is the angular veloc- 
ity (= ^pr)i and where we assumed -0 ~ 1 and a ~ 1 (a 
is viscous parameter) for the hard state. Clearly, for rea- 
sonable parameters, t pc and t comp are much smaller than 
t v i s , which means that inverse Compton scattering is an 
efficient cooling process in this situation. 

4-4- Qualitative Explanation for Anti- correlation be- 
tween Luminosity and Electron Temperature 

Next, we try to explain the anti-correlation between L 
and kT e quantitavely in the following simple discussion. 
The proton temperature (fcT p ) is assumed to be approxi- 
mately constant at ~ GA ^ ? " P , that is, the energy loss rate of 
protons is much smaller than viscous heating rate through 
the accretion. These protons will give their energy to elec- 
trons through collisions and the energy loss rate per unit 



volume is given as -^-r 
will get 



If 



is larger than Tp ~ , we 

rn~c. A 1 



tpe C*- 



(4) 



In a steady state, the heating rate from proton to elec- 
tron should balance with the cooling rate due to inverse 
Compton scattering, i.e. 



^nkTp 

tn& 



4kT. 



m c* 



(5) 



Using these equation (4) and (5), we can work out the anti- 
correlation ofkT e ocL~= quantitatively. This is attributed 
to the fact that the radiation mechanism in the hard state 
is due to inverse Compton scattering. 

Zdziarski (1998) suggested that the high energy cut- 
off depends on the luminosity from a theoretical point 
of view of thermal Comptonization by a more precise 
model. The relations kT c cx L~ 2 / 7 and kT c oc L -1 / 6 were 
derived for the advection dominated and cooling domi- 
nated cases, respectively. In comparison with this predic- 
tion, our result, kT e oc £~°- 24±0 06 , j s c l ose to these values. 

Furthermore, the maximum luminosity in the low/hard 

3 7 

state is predicted as L max ~ 0.15ysa5L Edd in the model 
of Zdziarski (1998). The highest luminosity observed in 
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2002, ~0.29LEdd) agrees within a factor of 2 with this 
value. 

In recent years the X-ray spectra of the hard state have 
also been discussed in the framework of jet outflows. In 
the most recent version of this model (Markoff, Nowak, 
& Wilms 2005) the hard X-ray spectra are dominated 
by synchrotron (at the low end) and synchrotron self- 
Compton (SSC) emission (at the high end) from the mag- 
netized plasma at the base of a jet, which shares many 
properties with other proposed types of Comptonizing 
coronae. In the simplest jet SSC interpretation, the 
observed decrease of the energy cut-off with luminosity 
means that either the size of the base of the jet is in- 
creasing or that the electron temperature is decreasing as 
the luminosity goes up. Changes are also expected in the 
synchrotron cut-off, which should decrease with luminos- 
ity, because of a higher cooling rate. There is likely an 
interplay between the synchrotron and the SSC compo- 
nents that needs to be worked out in more detail before 
the jet model can be applied to our data. 

4-5. Difference between bright hard state and common 
low/hard state 

Differences between the bright hard state (i.e. ^0.2 
L^dd) and the more common low/hard state have not been 
widely discussed in the literature. From an observational 
point of view, we found that the spectral properties in the 
bright hard state were not very different from the typical 
low/hard state except for one point: the high energy cut- 
off in the bright hard state is considerable lower (E cut «50 
keV) than in the low hard state (E cut «200 keV). The tim- 
ing properties are also suggested to be almost the same, 
with a moderate increase of the frequencies with lumi- 
nosity (Belloni et al. 2005). A clear correlation between 
the radio and X-ray flux has been established in the hard 
state of GX 339-4 (Corbel et al. 2003), but this relation 
has only been measured up to a luminosity of ~ 6 x 10 37 
erg s , i.e. just up to the luminosity above which we see 
a decrease in the energy cut-off. 

The low/hard state has been considered to correspond 
to an accretion-disk solution of the optically thin ADAF 
(Esin et al. 1997). Observed luminosity (~ 0.2LEdd) in 
the bright hard state can only be explained by this solu- 
tion (0.4a 2 LEdd) if we assume a large viscosity parameter 
a ~ 1. It is difficult to estimate an actual value of a, but 
recent three-dimensional MHD simulations suggest that 
a takes a value of reasonably 10~ 2 to 10 _1 (Hawley & 
Krolik 2001). Furthermore, considering the presence of 
hysteresis in the state transition, i.e. the fact that the 
bright hard state is only in the rise phase of the outburst, 
the bright hard state may suggest a presence of another 
solution like luminous hot accretion flow (LHAF: Yuan 
2001) rather than an ADAF. Further investigations in- 
cluding the radio/X-ray correlations are needed for the 
bright hard state through multi- wavelength observations. 

4-6. Constancy of the Photon Index Variations 

Another important result of our systematic study is that 
the slope of photon index remains almost constant at 1.4- 



1.7 regardless of the source intensity. This value is typical 
for the hard state of black hole candidates. From the point 
of view of thermal Comptonization, the photon index is 
correlated with Compton y parameter (= ^max(T,T 2 ); 
Sunyaev & Titartuck 1980). Thus, the obtained y pa- 
rameter distribution is remarkably close to unity for more 
than an order of magnitude as seen in Figure 8. This is 
consistent with power-law shape of the spectrum in the 
hard state. In the thermal Comptonization model, the y 
parameter should always be close to be unity when there 
is a large number of soft photons (Shapiro, Lightman, and 
Eardly 1976). The constancy of the y parameter means 
that the ratio of Comptonization luminosity to that in 
the seed soft photons is constant, indicating a constant 
geometry of the corona. 
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Fig. 8. Relation between the luminosity(L) and y parameter 
obtained by COMPST model. The different symbols corre- 
spond to the epoch defined in Figure 1. 



5. Summary 

We have analyzed 200 RXTE observations of the 
black hole candidate GX 339-4 in the hard state during 
its outbursts in 1996-2005. The observed luminosities 
ranged was from l.OxlO 37 erg s _1 to 2.1xl0 38 erg s _1 , 
assuming a source distance at 8 kpc. All the broadband 
spectra are well-explained by both cut-off power law 
(CPL) model and thermal Comptonization (COMPST) 
model with a spherical geometry. The photon index 
(in the cut-off power-law model) ranged from 1.4-1.7 
regardless of the luminosity. The Compton y parameter 
remained also close to unity regardless of the intensity 
variations, which is consistent with its power-law like 
spectral shape. The high energy cut-off in the CPL 
model ranges from 40 to 280 keV or more, and electron 
temperature in the COMPST model from 16 to 40 keV. 
Both parameters depend on the luminosity. We found a 
clear anti-correlation of E cut oc £~°- 70±0 - 06 between L 
and E cut or kT c oc £-0-24±o.06 between L and kT c when 
the L is 7xl0 37 erg s _1 . This former anti-correlation 
is quantitatively explained by the scenario that the 
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heating rate from protons to electrons balances with 
cooling rate of inverse Compton scattering. The cooling 
time scale of inverse Compton scattering obtained by 
spectral parameter is much smaller than viscous time 
scale, suggesting that the inverse Compton scattering is 
dominant in the radiation mechanisms. Further detailed 
studies of a high energy cut-off using INTEGRAL and 
Suzaku mission or future high energy missions will be 
valuable and interesting. 

The authors would like to thank the referee, Dr. Aya 
Kubota, for her valuable suggestions and comments. They 
also thank Sera Markoff for discussions on our results 
with regard to jet models for the hard state. This re- 
search has made use of data obtained through the High 
Energy Astrophysics Science Archive Research Center 
Online Service, provided by the NASA/Goddard Space 
Flight Center. 
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Table 2. Results of spectral fitting (I). 
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4.20+° 


41 
48 


1 


4±0 


1 


xlO" 


-9 


30. 91+." 


.20 
12 


4.21+° 


53 
76 


7 


1+0.3 
-"--0.2 


x 10~ 


9 


27.43+ 2 


15 
01 


4.61+° 


26 
27 


(9 


7±0.2 


xlO" 


-g 


24.83+° j 9 


4.76+.° 


09 
09 


1.7±0.04)xl0 


-8 




4.90+° 


12 
12 


1.9±0.1) 


xlO 


-8 


19.21+.°, 


76 
71 


5.09+.°, 


10 
15 


(2.1±0.1 


xlO" 


-8 


19 95+"- 85 


5.02±0 


40 



MAXt xVd.o.f 



40420 
40068 
40068 
20181 
20181 
20183 
20181 
20183 
20183 
20181 
20056 
20056 
20056 
20056 
20056 
20056 
20056 
20181 
20181 
20181 
20181 
20181 
20181 
20181 
20183 
20183 
20181 
20181 
20183 
20183 
20183 
40108 
40108 
40108 
40108 
40108 
40104 
60705 
70109 
70109 
70110 
70110 



01-01-00 
05-01-00 
05-02-00 
01-01-01 
01-01-00 
01-01-00 
01-02-00 
01-02-00 
01-02-01 
01-03-00 
01-01-00 
01-02-00 
01-03-00 
01-06-00 
01-06-01 
01-07-00 
01-08-00 
01-04-00 
01-05-00 
01-06-00 
01-06-01 
01-07-01 
01-07-02 
01-07-00 
01-03-00 
01-04-00 



1996-07- 
1996-10- 

1996- 10- 

1997- 02- 
1997-02- 
1997-02- 
1997-02- 
1997-02- 
1997-02- 
1997-02- 
1997-04- 
1997-04- 
1997-04- 
1997-04- 
1997-04- 
1997-04- 
1997-04- 
1997-05- 
1997-07- 
1997-08- 
1997-08- 
1997-09- 
1997-09- 
1997-09- 
1997-10- 
1997-10- 



01-08-000 1997-10 



01-08-00 
01-05-00 
01-06-00 
01-07-00 
01-03-00 
01-04-00 



1997-10 
1997-10 
1997-10 
1997-11 
1999-02 
1999-03 



01- 04-000 1999-03 

02- 01-00 1999-04 



02-02-00 
01-01-00 
01-55-00 
01-02-00 
01-01-00 
01-01-00 
01-02-00 



1999-04 
1999-04 
2002-04 
2002-04 
2002-04 
2002-04 
2002-04 



26 18:20-20:15 
17 02:12-05:07 
29 22:34-00:45 
03 15:56-19:09 
03 22:27-01:36 

08 14:20-20:48 
10 15:49-20:51 
44 00:09-06:43 
44 44:20-21:22 
17 18:29-00:12 
05 08:36-09:15 

10 11:47-12:28 

11 13:26-14:06 

17 23:26-00:01 

18 00:01-00:09 

19 22:20-23:21 

22 21:53-22:29 
29 09:23-13:51 
07 08:16-13:01 

23 04:23-07:29 
23 11:04-14:46 
19 22:09-22:37 

19 23:57-00:32 

20 01:30-05:58 
22 03:01-05:52 
25 03:22-06:16 
28 03:00-09:46 
28 09:46-11:46 
28 18:08-22:13 
31 20:35-23:48 
03 20:35-23:48 

12 21:40-01:54 
03 21:28-00:39 
03 15:07-21:28 
02 13:05-17:46 
22 10:46-17:09 

27 04:18-08:41 

02 44:24-44:54 

03 48:46-21:37 
07 13:07-14:07 

09 05:42-06:21 

10 10:11-10:22 



6.5±0.3 
1.7±0.1 

2.8 



o.i 

0.2 

6.6±0.3 
6.4±0.3 
6.1±0.2 
5.9±0.2 
5.9±0.3 
5.9±0.2 
5.8±0.3 
5.6±0.3 
6.2±0.3 
6.3±0.3 
6.5±0.3 
7.4 
6.5±0.3 
6.0±0.3 
3.9±0.2 
1.9±0.1 
4.9±0.3 
4.8±0.2 
5.7±0.3 
5.8±0.4 
5.8±0.3 
5.4±0.2 
4.6 



+0.3 
0.2 

4.1±0.2 

3 9+°' 3 

o.»_ 2 

4.0±0.2 
3.6±0.2 
3.6±0.2 
2.6±0.2 
3.2±0.2 
3.2±0.2 
3.7±0.2 
1.7±0.1 
1.4±0.1 
7.5±0.5 
L.l 



+ 0.04 
0.03 

1.8±0.1 
2.0±0.1 
2.2±0.1 



io-° 

XlO -9 
xlO -9 
XlO" 9 
xlO" 9 
XlO" 9 
x 10 _ 9 
x 10 — 9 
xlO" 9 
xlO- 9 
xlO -9 
xlO -9 
xlO- 9 
xlO" 9 
xlO -9 
XlO -9 
XlO" 9 
xlO" 9 
XlO" 9 
xlO- 9 
XlO -9 
XlO" 9 
xlO- 9 

io- 9 
io- 9 

9 

-9 
-9 
9 
,-9 



io- 

:10- 
:10~ 
:10~ 
:10~ 



x 

x 

X 
X 
X 
X 

x 

xlO- 9 
xlO -9 
xlO- 9 
xlO- 9 
xlO- 9 
xlO- 9 
x 10- 8 
xlO -8 
xlO- 8 
xlO -8 



167.7+. 24 2 
178.2 
243.4 



+ 31.1 
23.5 
+ 81.0 
-49.4 

168.6+ 2 |° 
162.5+ 3 5 8 
237.3+. 6 f 4 8 

, + 60.5 
-37.4 



179. 6Z 

165.7 



+ 47.5 
-31.2 



167. 5 1 



+ 55.0 
34.6 

> 198.4 
170 Q+ 601 
162.0+ 3 2 ° 
203.9+°°;° 

>232.7 
235.31^ 
241.2±|f| 
179.8+. 86 .; 2 
237.3+5, 9B 4 7 

, + 39.0 



27.6 
+ 60.4 

39.3 
+ 71.1 



215. T 



176.2; 
213.3 
204.3_ 4 
217.4+.63;! 

, + 85.9 
-48.4 

>210.7 
>254.3 

249 g+ 144 -° 
>225.1 
>223.0 
>254.9 
>226.8 
>174.6 

203.9+™° 
114.9+™ 
84.9+° 3 



75.47; 
79.9 



+ 7.8 
6.7 



1.52±0.03 

1 49+ 04 

1 50+°° 3 
1 - OU -0.04 

1.52±0.02 
1.52±0.03 
1.50±0.02 
1.49±0.02 
1.52±0.03 
1.47±0.02 
1.47±0.04 
1.54±0.04 
1.55±0.03 
1.51±0.03 
1.55±0.03 
1 59+ 04 

±-tJJ -0.06 

1.55±0.03 
1.56±0.04 
1.49±0.03 
1 49+ 03 

-"-■^ -0.04 

1.50±0.04 
1.50±0.03 
1.52±0.04 
1.54±0.04 
1.50±0.03 
1.55±0.03 
1.52±0.04 
1.50±0.03 
1.51±0.04 
1.48±0.04 
1.52±0.04 
1.51±0.03 
1.62±0.04 
1 54+ 02 

- L - J *-0.04 

1.54±0.04 
1 49+0 03 

-0.04 
, ,-,+0.04 
1 - ol -0.05 

1.52±0.05 
1.44±0.04 
1.36±0.02 
1.41±0.02 
1.44±0.03 
1.44±0.03 



1.45±0.24 

96+ ' 39 

u ' 3u -0.36 

96+ 34 
1.41±0.23 
1.46±0.25 
1.71±0.21 
1.61±0.22 
1.48±0.26 
1.85±0.21 
1.80+ 



-, + 0.32 
-0.31 

1.20±0.30 
1 3n+ ' 28 
1.43±0.28 
1.43±0.28 

1 29+ ' 45 

1 ' za -0.42 

1 54+ 028 
1.54±0.29 
1.41±0.29 

1 00+ 034 

0.31 

1.62±0.35 
1.37±0.27 
1.37±0.32 

1 qc + 0.35 
i " 3D -0.34 

1.87±0.28 
1.43±0.24 
1.39±0.30 
1.39±0.27 
1.11±0.35 
1.43±0.31 
1.06±0.30 
1.12±0.27 
1.56±0.29 
1.43±0.32 
1.23±0.33 
1.52±0.33 
1.131 
0.93 4 



, + 0.43 
-0.20 
3 + 0.41 
-0.40 

1.28±0.35 
1 94+ 019 
1 93+ 23 
1.80±0.25 
1.79±0.27 



80.8/79 
74.9/79 
64.3/79 
79.2/79 
62.0/79 
82.7/79 
63.1/79 
95.9/79 

115.6/79 
77.4/79 
64.8/79 
68.2/79 
79.0/79 
56.1/79 
68.2/79 

52.77/79 
71.6/79 
78.1/79 
65.8/79 
72.5/79 
53.0/79 
46.7/79 
58.4/79 
77.5/79 
58.6/79 
60.8/79 
81.3/79 
74.1/79 
73.3/79 
49.8/79 
56.1/79 
60.2/79 
64.1/79 
62.7/79 
52.0/73 
61.0/73 
57.7/73 
80.8/73 
79.9/73 
66.5/73 
62.9/73 
81.5/73 



0.71±0.19 
0.60±0.33 
55+ ' 28 
0.65±0.18 
0.69±0.20 
0.75±0.17 
0.73±0.18 
0.72±0.2iq 
0.79±0.JE 
0.60±0.g3 

o.73i°;| 

0.72±0.p 
0.83^10.^ 
1.03±0.^ 
0.84±0.^ 
1.00±0. 




i.06±o.;._ 

0.86±0.|f 
0.80±0.2) 
0.70±0.2 J i 
0.60±0.gl 
0.73±0.^) 
0.81±0.fe 
0.62±0.|£ 
0.65±0.^ 
1.13±0.27 
96+ 26 

U.»U_ Q 25 

0.60±0.25 
0.77±0.25 

n Sfi+0-3^ 

-0.33 

0.45 ' 
0.41^ 
0.26 4 
0.46 4 



-+0.35 
-0.34 
I +0.28 
-0.27 
; + 0.15 
-0.15 
;+0.19 
-0.19 

0.54±0.21 

-, + 0.23 
-0.22 



0.60 1 



90.1/79 
77.4/79 
61.6/79 
91.5/79 
70.7/79 
79.4/79 
73.3/79 
92.8/79 
111.1/79 
80.8/79 
68.2/79 
77.9/79 
83.1/79 
55.6/79 
73.3/79 
60.7/79 
73.5/79 
92.1/79 
80.5/79 
76.2/79 
55.2/79 
42.0/79 
47.3/79 
82.3/79 
62.8/79 
71.2/79 
74.8/79 
63.4/79 
58.0/79 
45.5/79 
51.9/79 
52.7/79 
58.0/79 
59.9/79 
66.9/73 
66.0/73 
69.8/73 
93.7/73 
140.9/73 
139.5/73 
111.5/73 
125.7/73 



Err. 



ors arc . 



%coiili(lciicc.Thc unabsorbed Huxis calculated for 2-200 kcV. 



Table 3. Results of spectral fitting (II). 



cut-off power-law 



(CPL) 



thermal 



Flux (crg/scc/cm^) 



Comptoiiization (COMPST) 



Observation ID 


Observation Time 


Flux (erg/ 


sec / cm 


70109 


01 


03 


00 


2002 


04 


13 


02 


06 


02 


30 


(2.3±0.1 


X10 


-8 


70110 


01 


03 


00 


2002 


04 


15 


03 


13 


03 


40 


(2.5±0.1 


xlO 


-8 


70110 


01 


04 


00 


2002 


04 


17 


02 


43 


03 


09 


(2.6±0.1 


xlO 


-8 


40031 


03 


01 


00 


2002 


04 


18 


02 


28 


03 


12 


(2.6±0.1 


xlO 


-8 


40031 


03 


02 


00 


2002 


04 


19 


02 


12 


02 


40 


(2.7±0.1 


xlO 


-8 


40031 


03 


02 


01 


2002 


04 


20 


02 


55 


04 


09 


(2.8±0.1 


xlO 


-8 


40031 


03 


02 


02 


2002 


04 


20 


22 


24 


00 


37 


(2.7±0.1 


xlO 


-8 


70110 


01 


05 


00 


2002 


04 


21 


11 


09 


11 


24 


(2.7±0.1 


xlO 


-8 


40031 


03 


02 


03 


2002 


04 


21 


22 


07 


03 


07 


(2.7±0.1 


xlO 


-8 


40031 


03 


02 


04 


2002 


04 


22 


21 


51 


23 


07 


(2.8±0.1 


xlO 


-8 


70110 


01 


06 


00 


2002 


04 


2:5 


11 


39 


11 


54 


(2.8±0.1 


xlO 


-8 


70109 


01 


04 


00 


2002 


04 


23 


13 


20 


18 


13 


(2.7±0.1 


xlO 


-8 


40031 


03 


02 


05 


2002 


04 


24 


01 


19 


01 


48 


(2.8±0.1 


xlO 


-8 


40031 


03 


02 


06 


2002 


04 


25 


01 


33 


01 


59 


(2.7±0.2 


xlO 


-8 


70110 


01 


07 


00 


2002 


04 


26 


21 


43 


22 


20 


(2.7±0.1 


xlO 


-8 


70109-01-05-01G 


2002 


04 


27 


07 


39 


08 


50 


(2.8±0.1 


xlO 


-8 


70109 


01 


05 


00 


2002 


04 


27 


09 


12 


10 


36 


(2.7±0.1 


xlO 


-8 


70109 


01 


05 


02 


2002 


04 


27 


10 


48 


11 


51 


(2.8±0.1 


xlO 


-8 


70110 


01 


08 


00 


2002 


04 


30 


10 


28 


11 


02 


(2.8±0.1 


xlO 


-8 


70110- 


01 


01 


10 


2003 


04 


11 


00 


19 


00 


36 


(2.0±0.2 


xlO 


-9 


80116 


02 


02 


00 


2003 


04 


12 


16 


45 


19 


56 


(2.0±0.1 


xlO 


-9 


80116 


02 


02 


01 


2003 


04 


13 


00 


55 


01 


3,2 


v —0.2 


xlO 


-9 


70128 


02 


03 


Of) 


2003 


04 


13 


05 


39 


11 


00 


(1.9±0.1 


xlO 


-9 


70128 


02 


03 


01 


2003 


04 


13 


14 


40 


19 


35 


(2.0±0.1 


xlO 


-9 


70110 


01 


02 


10 


2003 


04 


14 


03 


48 


04 


02 


(1.9±0.2 


xlO 


-9 


80116-02-02-02G 


2003 


04 


14 


06 


22 


09 


28 


(1.9±0.1 


xlO 


-9 


60705 


01 


61 


00 


2003 


04 


17 


23 


05 


00 


01 


(1.3±0.1 


xlO 


-9 


70110- 


01 


03 


10 


2003 


04 


18 


09 


39 


10 


06 


(1.4±0.2 


xlO 


-9 


80132 


01 


05 


00 


2004 


02 


21 


09 


16 


10 


10 


(1.1+0.1 


xlO 


-9 


80132 


01 


05 


01 


2004 


02 


21 


13 


10 


14 


17 


1.8+2-? 

— U. 1 


xio~ 


9 


80132 


01 


06 


00 


2004 


02 


22 


05 


34 


06 


29 


(2.0±0.2 


xlO 


-9 


80132 


01 


06 


01 


2004 


02 


22 


20 


51 


22 


14 


(1.9±0.1 


xlO 


-9 


80132 


01 


07 


01 


2004 


02 


23 


20 


30 


21 


52 


9 9 + 0-1 
^-0.2 


XlO" 


9 


90109 


01 


01 


00 


2004 


02 


23 


23 


43 


04 


11 


(2.1±0.1 


xlO 


-9 


80102 


04 


58 


01 


2004 


02 


24 


14 


18 


14 


46 


(2.2±0.2 


xlO 


-9 


80132 


01 


07 


00 


2004 


02 


26 


17 


56 


20 


32 


(2.5±0.2 


xlO 


-9 


80102 


04 


59 


00 


2004 


02 


27 


10 


37 


11 


15 


(2.9±0.3 


xlO 


-9 


80132 


01 


08 


00 


2004 


02 


27 


19 


08 


20 


42 


(2.6±0.2 


xlO 


-9 



E out (keV) 

63.3+|* 
60.9^' 9 7 

+ 2.9 
-2.7 
+ 2.8 
-2.6 

58.8^21 
58.1+^ 
60.8+a'i 
60. 9j" " 
60. 1 4 



photon index MAXr x 





02 


2.1012 


21 


58 


4/73 


(2.2 + 0.1 


xlO 


-8 


1.44±0 


02 


1.83+0 


21 


46 


8/73 


(2.4+0.1 


xlO 


-8 


1.45+0 


02 


2.06+0 


20 


80 


3/73 


(2.5+0.1 


xlO 


-8 


1.46+0 


1)2 


1.77+0 


20 


53 


2/73 


(2.5+0.1 


x 10 


-8 


1.47+0 


02 


1.75+0 


19 


62 


1/73 


n e+0.04 
"■° — 05 


x 10" 


-8 


1.45+0 


02 


1.97+0 


If) 


69 


4/73 


9 R +0.'04 
Z - D _0 01 


x 10" 


-8 


1.44+0 


1)2 


2.1512 


16 


93 


1/73 


9 -+0.05 
— 0.04 


x 10" 


-8 


1.46+0 


02 


1.86+0 


22 


65 


1/73 


(2.6+0.1 


xlO 


-8 


1.47+0 


01 


1.96+0 


13 


02 


9/73 


(2.6+0.4 


xlO 


-8 


1.47+0 


1)2 


1.91+0 


17 


74 


9/73 


9 fi +0.05 
Z - D — 0.04 


x 10" 


-8 


1.50+0 


03 


1.67+0 


27 


62 


4/73 


(2.7+0.1 


xlO 


-8 


1.4712 


01 

01 


1.9712 


12 
12 


127.7/73 


9 ,5+0.03 
Z - D _0.04 


x 10" 


-8 


1.48+0 


02 


1.83+0 


20 


68 


0/73 


(2.6+0.1 


xlO 


-8 


1.46+0 


04 


1.85+0 


36 


60 


6/73 


(2.7+0.1 


x 10 


-8 


1.46+0 


03 


1.90+0 


27 


55 


4/73 


(2.7+0.1 


xlO 


-8 


1.50+0 


02 


1.72+0 


18 


47 


3/73 


(2.6+0.1 


xlO 


-8 


1.50+0 


1)2 


1.73+2 


18 
17 


71 


1/73 


(2.6+0.1 


xlO 


-8 


1.50+0 


1)2 


1.7112 


17 
17 


77 


1/73 


(2.6+0.1 


xlO 


-8 


1.53+0 


02 


1.72+0 


19 


78 


0/73 


(2.6+0.1 


xlO 


-8 


1.66+0 


08 


1.3112 


63 


65 


6/74 


i.9l2;? 


K 10" 


9 


1.66t° 


1 
05 


1.37+0 


38 


56 


9/74 


(1.9+0.1 


xlO 


-9 


1.631° 


03 


1.5112 


10 


65 


2/74 


(2.1+0.1 


xlO 


-9 


1.59+0 


03 


1.49+0 


29 


68 


6/74 


(1.8+0.1 


x 10 


-9 


1.59+0 


04 


1.54+0 


36 


57 


3/74 


(1.9+0.1 


xlO 


-9 


1.64+2 


02 
08 


1.4012 


60 
54 


59 


9/74 


(2.1+0.2 


x 10 


-9 


1.59+0 


04 


1.39+0 


35 


04 


4/74 


(1.8+0.1 


x 10 


-9 


1.55+0 


06 


1.23+0 


48 


60 


0/74 


(1.3+0.1 


xlO 


-9 


1.48+2 


10 


2.29+0 


79 


48 


7/74 


(1.3+0.1 


x 10 


-9 


1.46+2 


05 
07 


1.1812 


53 
50 


73 


7/74 


(1.8+0.1 


xlO 


-9 


1.4312 


05 
06 


1.27+0 


48 


70 


2/74 


(1.7+0.1 


xlO 


-9 


1.4612 


1 
06 


0.89+2 
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Table 4. Results of spectral fitting (III). 



cut-off power-law (CPL) thermal Comptonization (COMPST) 
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E cM (kcV) 


photon index 


MAX 


r 


\ J 


/d.o.f 


Flux (erf 


,-/ 


sec / cm") 


kT (keV) 


T 




MAXt x'Vd-O-f 


80132 


01 


08 


01 


2004 


02 


28 


07 


51 


08 


53 


(2 
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Table 5. Results of spectral fitting (IV). 



cut-off power-law (CPL) thermal Comptouizalion (COMPST) 



Observation ID 


Observation Time 


Flux (crg/scc/cm^) 
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Table 6. Results of spectral fitting (V). 



cut-off power-law (CPL) thermal Comptonization (COMPST) 



Observation ID 
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